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Abstract: We analyse the sensitivity of geosynchrotron radio emission from inclined extensive air show- 
ers to the energy and mass of primary cosmic rays. We demonstrate that radio emission measurements 
at suitable lateral distances can infer both the number of electrons and positrons in the shower maximum 
and the atmospheric depth of the maximum on a shower-to-shower basis. Alternatively, measurements at 
a fixed lateral distance but in two different observing frequency bands yield comparable information. An 
RMS error of 5% in the determination of the number of electrons and positrons at shower maximum can 
be achieved. Through the determination of these quantities, geosynchrotron radiation provides access to 
the energy and mass of primary cosmic rays on a shower-to-shower basis. 
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Introduction 

It is our current understanding that radio emis- 
sion from extensive air showers is dominated by 
"geosynchrotron radiation" [3] emitted by sec- 
ondary shower electrons and positrons being de- 
flected in the earth's magnetic field. The geosyn- 
chrotron model has by now been implemented in a 
sophisticated Monte Carlo code called REAS2 [5, 
6] , which itself relies on CORSIKA [ 1 ] for the sim- 
ulation of the relevant extensive air shower (EAS) 
properties. In this article, we analyse the sensi- 
tivity of REAS2-simulated geosynchrotron radio 
emission on the primary cosmic ray energy and 
mass. The determination of these parameters on 
a shower-to-shower basis is one of the main goals 
of measuring radio emission from EAS, and we 
demonstrate that radio measurements of inclined 
showers indeed provide relatively direct access to 
these parameters. 

Methodology 

Results gathered so far point to inclined air show- 
ers with zenith angles above 45° as particularly 
promising targets for radio measurements of EAS 
[7]. Inclined air showers have a large radio "foot- 
print" [4] and thus allow radio antennas to be 
spaced relatively far apart, an important prerequi- 



site for the instrumentation of large effective ar- 
eas at moderate cost. Concentrating on the en- 
ergy range relevant to the Pierre Auger Observa- 
tory, we have thus performed an analysis of air 
showers with 60° and 45° zenith angles at pri- 
mary particle energies of 10 18 eV, 10 19 eV and 
10 20 eV. For each of these zenith angles and en- 
ergies, we have simulated 25 iron-induced and 25 
proton-induced air showers. 25 gamma-induced 
air showers per zenith angle were simulated for 
10 18 eV and 10 19 eV, but not for 10 20 eV, where 
pre-showering in the geomagnetic field would have 
to be taken into account. The simulation chain 
consisted of a CORSIKA 6.502 run using the 
QGSJETII-03 and UrQMD 1.3.1 interaction mod- 
els, Argentinian magnetic field, fixed azimuth an- 
gle (showers coming from south), optimised 
thinning, and an observer height 1400 m a.s.l. fol- 
lowed by a REAS2 simulation with antenna loca- 
tions between 25 m and 925 m from the shower 
core (in ground-based coordinates); cf. also [5]. 
For each of the simulated radio events we then de- 
termined the peak amplitude of the electric field 
pulses filtered using idealised rectangle filters from 
16 to 32 MHz, 32 to 64 MHz and 64 to 128 MHz, 
respectively. In the following, we discuss the 60° 
zenith angle case. The qualitative behaviour at 45° 
is completely analogue. Other shower azimuth an- 
gles do not change the qualitative behaviour either. 
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Figure 1: Peak radio field strengths of 10 19 eV 
showers with 60° zenith angle as measured in the 
flat region 275 m north of the shower core. 
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Figure 2: Peak radio field strengths of 10 19 eV 
showers with 60° zenith angle as measured in the 
steep region 725 m north of the shower core. 



Signal information content 

The lateral slope of the radio signal is known to 
exhibit a dependence on the depth of the shower 
maximum (X max ) and consequently, contains in- 
formation on the primary particle energy and mass 
[4, 2]. To date, however, there had been no de- 
tailed investigation how this information content 
could be exploited in practice. 
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Figure 3: The ratio of field strengths in the steep 
region vs. the^?af region of 10 19 eV showers with 
60° distinguishes between primary particle types. 

The detailed analysis presented here reveals that 
for a given shower zenith angle, a suitable lateral 
distance exists where the radio signal is relatively 
independent of JT max (hereafter called flat region), 



in analogy to the surface detector quantity 5(1000) 
of the Pierre Auger Observatory. This behaviour is 
illustrated for the 60° zenith angle case in Fig. 1 : 
the electric field strength at 275 m north from the 
shower core is relatively constant regardless of pri- 
mary particle type and shower -X" max . In contrast, 
at a distance of 725 m north as shown in Fig. 2, 
there is a clear dependence of signal strength on 
X max (hereafter called steep region). 

A combination of measurements in these two re- 
gions can therefore differentiate between different 
types of primaries (Fig. 3). 

Signal scaling with iV max 

When the measured radio field strength is divided 
by the number of electrons plus positrons in X max , 
hereafter called N max , it becomes clear that the 
radio signal scales linearly with iV max . The rea- 
son for this is that most of the radio emission 
stems from the particles close to the shower max- 
imum [5]. (The intensity of optical fluorescence 
light, in contrast, scales with the calorimetric en- 
ergy deposited in the atmosphere by the shower.) 
The clean linear scaling is illustrated in Fig. 4, 
where the results from all energies and particle 
species yield an approximately constant electric 
field strength per iV max in the flat region. The elec- 
tric field strength per iV max in the steep region is 
also constant over the different particle types and 
energies for a given X max . 
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Figure 4: In the flat region at 275 m north, the 
radio pulse height per iV max over all energies and 
particle types is approximately constant. 




600 650 700 750 800 850 900 950 100010501100 
Xmax [g cm" 2 ] 

Figure 5: In the steep region at 725 m north, the 
radio pulse height per iV max over all energies and 
particle types is constant for a given X max . 
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Figure 6: A measurement in the flat region at 
275 m north directly yields the shower N max . The 
line denotes a linear functional dependence. 
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Figure 7: The ratio of peak field strengths in the 
flat region and steep region yields direct informa- 
tion on the shower X max . 



Optimum parameter determination 

Histogramming the distribution of electric field 
strengths per iV max over all particle species and 
energies yields an RMS of only 5% in the flat re- 
gion (Fig. 8). In spite of shower to shower fluc- 
tuations, radio measurements could thus determine 
the N max of individual EAS to very high precision. 
The position of the flat region for a given zenith 
angle is determined by the observing frequency 
band. As illustrated in Fig. 9, for 60° zenith angle 
and 32 to 64 MHz observing frequency, N mlLX can 
be determined to highest precision around 275 m 
distance to the north (or south). At 16 to 32 
MHz, 5% precision in the A^ nax determination can 
be reached anywhere between around 200 m and 
500 m from the core. Low frequencies are thus 



particularly well suited for an energy determina- 
tion of EAS primaries, if the technical difficulties 
involved with measurements at these frequencies 
can be overcome. In contrast, at 64 to 128 MHz 
the range for iV max determination becomes much 
smaller and resides at smaller distances. (The pe- 
culiar steps in the 64 to 128 MHz RMS curve point 
to coherence effects starting to play a significant 
role at these high frequencies.) For the determina- 
tion of X max using measurements in the steep re- 
gion, a lateral distance should be chosen that shows 
a large RMS spread in electric field strength per 
Nmax- For 32 to 64 MHz and 60° zenith angle, 
725 m constitutes a suitable compromise between a 
good handle on the X max value and detectable ab- 
solute signal levels. (The absolute signal strength 
drops quickly with lateral distance, cf. [4]). 
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Figure 8: For 60° zenith angle showers, the RMS 
spread in peak field strength per N max in the flat 
region over all energies and particle types amounts 
to ~ 5%. In principle, N max can thus be inferred 
with 5% uncertainty. 
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Figure 9: Determination of the RMS uncertainty in 
peak field strength per N max as a function of lat- 
eral (ground-coordinate) distance from the shower 
core. The flat region shifts to larger lateral dis- 
tances at lower frequencies. 



Conclusions 

Simulations of geosynchrotron radio emission with 
the REAS2 Monte Carlo code reveal that informa- 
tion contained in the lateral profile of the radio 
signal can be exploited for a direct determination 
of iV max and X max on a shower-to-shower basis. 
These quantities can in turn be related to the en- 
ergy and mass of the primary cosmic ray particle. 

For a given zenith angle and observing frequency 
band, distinct distance regimes denoted here as the 
flat region and the steep region can be identified. 
Measurements in the flat region directly yield the 
JV max of the shower, with an RMS deviation of 
only 5%. A comparison of peak field strengths 
in the steep region and flat region for a given ob- 
serving bandwidth provides a direct estimate of the 
shower X max . Alternatively, measurements at a 
fixed distance but in two different observing fre- 
quency bands provide comparable information. 

The position of the flat region in case of 60° zenith 
angles for the observing frequency band from 32 
to 64 MHz lies at approximately 275 m and would 
thus require an antenna spacing of that order. (The 
given number is valid for antenna positions along 
the shower axis; in the perpendicular direction the 
scales are smaller due to azimuthal asymmetries of 
the radio footprint.) If technical problems arising 
in measurements at lower frequencies such as 16 
to 32 MHz can be overcome, antenna spacings of 



up to 500 m will allow measurements in the flat re- 
gion. The qualitative behaviour at 45° zenith angle 
is analogue, yet at lower lateral distances. At larger 
zenith angles, the scales will be larger, but con- 
sistent simulations require the implementation of 
curved atmospheres for zenith angles much larger 
than 60°, which is currently being prepared. 
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